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Magnetic ions such as Fe, Mn and Co with locaUzed spins may be adsorbed on the surface of 
topological insulators like Bi2Se3. They form scattering centers for the helical surface states which 
have a Dirac cone dispersion as long as the local spins are disordered. However, the local density 
of states (LDOS) may be severely modified by the formation of bound states. Commonly only 
elastic scattering due to normal and exchange potentials of the adatom is assumed. Magnetization 
measurements show, however, that considerable magnetic single ion anisotropics exist which lead to 
a splitting of the local impurity spin states resulting in a singlet ground state. Therefore inelastic 
scattering processes of helical Dirac electrons become possible as described by a dynamical local self 
energy of second order in the exchange interaction. The self energy influences bound state formation 
and leads to significant new anomalies in the LDOS at low energies and low temperatures which we 
calculate within T-matrix approach. We propose that they may be used for spectroscopy of local 
impurity spin states by appropriate tuning of chemical potential and magnetic field. 

PACS numbers: 73.20.-r, 73.50.-h, 73.50.Bk, 75.30.Hx 



I. INTRODUCTION 

The warped Dirac cone dispersion and helical spin po- 
larization of surface states in topological insulators (TI) is 
by now well understood, in particular from ARPES inves- 
tigations [11, i2j . These states are topologically protected 
as long as time reversal invariance is preserved The 
latter may be broken by an applied magnetic field or by 
adsorption of magnetic adatoms in sufficient concentra- 
tion such that they exhibit long range order In this 
case the Dirac electrons acquire a mass gap proportional 
to the size of the magnetic field or magnetization. 

Even without long range spin order magnetic adatoms 
still have a profound efi^ect on the surface states. While 
there is no gap and at sufficient distance from the im- 
purity the unperturbed Dirac spectrum is recovered, the 
local density of states (LDOS) may be severely modi- 
fied due to formation of low energy bound states at the 
impurity site [6,]. In addition the spatial dependence of 
the DOS away from the impurity may be Fourier trans- 
formed to give the quasiparticle interference (QPI) spec- 
trum 0411 ■ It is dominated by the special wave vectors 
that connect nested Fermi surface pieces not related by 
backscattering. Therefore the QPI contains considerable 
information on Fermi surface geometry that may be com- 
pared to ARPES results [13, El- 

In all cases it is assumed that the scattering mecha- 
nism is elastic and no internal degree of freedom may 
be excited at the impurity site by the scattering. In 
the case of magnetic impurities it is commonly assumed 
that the moment is polarized along a fixed direction 
(the ensuing gap opening due the associated time 
reversal symmetry breaking is neglected). However 
this is an oversimplified picture of the real situation. 
Firstly in zero field and without order the spins are not 
polarized but fluctuate. Secondly due to the symmetry 
breaking the directions parallel and perpendicular to 
the surface plane are not equivalent and therefore 



single-ion anisotropics for the magnetic moment may 
be introduced. There is direct experimental evidence 
for this effect coming from rnagnetization measurements 
on Fe adsorbed on Bi2Se3 which suggest that the 

in-plane magnetization is considerably larger than the 
out-of-plane magnetization. This means that there is 
an easy-plane spin anisotropy which splits the quantum 
states of the the S" = 2 Fe impurity spin into three 
levels. Such anisotropy was also observed in the case 
of Fe adatoms on graphene d3i|. Therefore in second 
order of the exchange constant there will be inelastic 
scattering processes of helical Dirac electrons from those 
local impurity levels. 

The internal impurity excitation degrees of freedom 
can profoundly modify the low energy scattering matrix 
and therefore the bound state formation and LDOS. 
Since the excitation energies of the impurity depend 
on the field this should also lead to a field dependence 
of the LDOS even in the case when the ground state 
of the magnetic adatom is a singlet. Furthermore 
the LDOS will be temperature dependent due to the 
thermal occupation of the impurity levels. One may 
also expect that these effects depend on the distance of 
the chemical potential from the Dirac point as long as 
it is comparable to the splitting energy. We will study 
these dynamical impurity self energy effects in this work 
and show that they open interesting perspectives for 
the STM experiments on topological insulators with 
magnetic adatoms. 

We mention that the effect of inelastic scattering 
on LDOS properties has been investigated before in a 
completely different context: In superconductors bound 
states due to impurity scattering may appear within 
the superconducting gap The physics of these 

bound states is strongly modified when the impurities 
are dynamic, i.e, exhibit a local phonon mode [l5| or in- 
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FIG. 1. Left: Energy level scheme of localized spin states split 
by the anisotropy term in Eq. ([T]) with A — Ko and A' — 4Ko. 
Right: Second order self energy for Dirac electrons (full line) 
due to scattering from local spin states a, (3 (broken line). 



tramolecular excitations [l6| . This problem has a certain 
formal analogy to the one studied here for a completely 
different and non-superconducting system. Spin-inleastic 
excitations and their influence on Friedel oscillations for 
normal metals were investigated in Refd 

In Sec. HI] we will first discuss the magnetic anisotropy 
and its connection to quantum states of adatoms with 
spin. Then we introduce the model for normal and mag- 
netic impurity scattering in Sec. IIIII The scattering ma- 
trix and dynamical self energy for the coupled problem 
are treated in Sec. IIVI The zero temperature self energy 
will be calculated analytically in Sec. |V] for comparison 
with numerical results and the resulting LDOS for vari- 
ous physical regimes will be discussed in Sec. I VII Finally 
Sec. IVIII gives the summary and conclusion. 



II. MAGNETIC STATES OF Fe SURFACE 
ADATOMS 

Various experiments have tested the robustness of TI 
surface states against doping with magnetic ions. Bulk 
doping Bi2Se3 with 1% Mn is sufficient to open a mag- 
netic mass gap 0- Surface deposition of Fe on the (111) 
selenium surface plane is also reported [l^ to lead to mas- 
sive Dirac dispersion of surface states. On the other hand 
these states are insensitive under sufficiently low doping 
of the same surface with magnetic Gd atoms [l^. Ap- 
parently sufficient surface doping is required to produce a 
static effective field that leads to the breaking of time re- 
versal symmetry necessary for creating a gap. In addition 
there is also a more local effect on the electronic structure 
close to the dopant atom which comes from repeated scat- 
tering of Dirac particles at the dopant site and may lead 
to resonant enhancement of the LDOS at low energies. 
To describe this effect it is necessary to have a valid model 
for the possible spin states of the adatom. Due to the spa- 
tial symmetry breaking at the surface these states may be 
split into a sequence of multiplets and lead to anisotropic 
magnetic properties at the surface. The latter have been 
reported for Fe doped (111) surfaces of Bi2Se3. From 
magnetization measurements an in-plane anisotropy was 



found [iJl . The adsorbed Fe atom is in a S* = 2 high-spin 
state with moment m = OsSfj-B = 4/iB and the orbital 
moment nearly quenched [12l | . The anisotropic single ion 
spin Hamiltonian is given by 



(1) 



Here — {K/ S'^) is the magnetic anisotropy constant. 
It may be directly determined from the full magnetiza- 
tion curves. It also determines the low-field susceptibility 
anisotropy. The reduced susceptibility x = xligst^sY 
for a field H forming an angle 9 with the surface normal 
z is given by 
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The dimensionless anisotropy constant k is then deter- 
mined by K = X||/x_L — 1 where || , _L refer to = ^,0 field 
direction, respectively. In absolute units the anisotropy 
was determined as Kq — 0.4 meV from the magnetiza- 
tion curves jl^l • Due to the positive sign the Fe moment 
is preferentially aligned with the surface plane. Microc- 
socpically it splits the five \Sz) local spin (S = 2) states 
into a nonmagnetic singlet ground state and two excited 
magnetic doublets \Sz) = |0), | ± 1), | ± 2) with energies 
Ea=0,A = Ko, A' = 4:Ko (a = 0, ±1, ±2) respectively 
(Fig. [T]). To treat the scattering of Dirac particles from 
the impurity spin states we describe the latter by a con- 
venient pseudo fermion representation [20l . [2l| : 
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(3) 



where the pseudofermion operators al^{i) create the five 
Fe spin states at impurity site i. Here Q is a Lagrange 
parameter used to project out unphysical states as ex- 
plained in Appendix |B] They correspond to the occu- 
pation of two or more different local spin states at the 
same site. For Fe on Bi2Se3 the total splitting is 1.6 
meV which is quite small compared to the surface band 
cutoff Ec — 0.25 eV. In our model calculations we will 
assume that the splitting may be considerably larger of 
the order A/Ec — 0.1 to exhibit clearly the qualitative 
effect of inelastic scattering on the LDOS. Such splitting 
energies seem entirely possible with magnetic 4f adatoms 
on Bi2Se3. 



III. MODEL FOR NORMAL AND MAGNETIC 
IMPURITY SCATTERING 

In this work we will be mostly interested in the ef- 
fects of inelastic magnetic scattering on the helical sur- 
face state LDOS at relatively low energies comparable to 
the splitting energy A but much smaller than the bulk 
gap. Therefore these states may be described by the 
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FIG. 2. (Color online) (a) real part 
of conduction electron Green's function 
Go (a;)' from analytical (full line) and nu- 
merical (broken line) calculation. The 
latter is performed with a soft cutoff of 
width j/Ec = 0.05 at E^._ (b) Corre- 
sponding imaginary parts Gq{uj)" . The 
DOS of Dirac half cones (r = ±1) is given 
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isotropic Dirac cone Hamiltonian neglecting the warping 
effect which is negligible in this region. It is given by 

^c = E*k[<kx<T).z]vI'k = ^vl-[/ik^j^ (4) 

k k 

with ^E*]^ = (cj^-f-, cjj.^) denoting the surface states. It is 
diagonalized by the helical eigenstates ^j^. = (7^_,_ , 7k- ) 
where r = ±1 is the helicity. Explicitly /i^ = EyT^ with 
the helical eigenvalues £'kr — tE\^ given by = u|k| 
where v is the velocity of Dirac particles. Here r — ±1 
denote the upper and lower Dirac cone respectively. The 
eigenstates are obtained by 

*k = P^*k; Pk=(_J-.,, '"7*") (5) 

where 0k = tan^^(fcj,/fca:) is the azimuthal angle of k. 

The scattering of these surface states from normal and 
magnetic impurities is described by contact interactions 
of strength Uq and Jq at the impurity sites R^: 

[/(r-R,) = [/o(5(r-R,) 

F(r-R,) = ijo<T.5(5(r-R,) (6) 

Their Fourier transforms are momentum independent 
with ?7kk' = Uo and Vkk' = jJof ■ S. The total scat- 
tering Hamiltonian Hj = Hjj + Hj is then given by 

Ha-^oEE4+q.Ck.'e^'»^^ (7) 

ikq a(5{T{7' 

The matrix elements result from the pseudo-fermion rep- 
resentation for spin operators which may be expressed 

as 

= ^{a\S^\l3)a)^{i)ai3{i) 

M^^p ^ {a\S,\l3) ■ (T,,, (8) 

Here a, /3 denote the split impurity spin states. The total 
model Hamiltonian is then given by "H Ha + Hd + 
Hi- It is necessary to transform H to the helical basis 
states created by 7^^ using the unitary transformation 



in Eq. ([S]). We separate the mean field term ^ (S) of 
the exchange part and add it to the potential scattering 
Then we obtain 

'^f^ = E E k + q)7,l+q.7k.'e'''^- (9) 

i kqr 

i kqrr' 

For the mean field part we have, denoting k' = k + q 

C/^^'(kk') = UoU^^'ikk') + i7ii7r'(kk') (10) 

The explicit form of matrix elements C/^^ (kk') and 
A/^^ (kk') is given in Appendix [A] Assuming that the 
mean field spin polarization is oriented along z the corre- 
sponding scattering strength is given by f/i = {Jo/2){Sz)- 
The spin polarization can appear as result of spontaneous 
FM order of impurity spins for sufficient surface cover- 
age or through an applied field H. We note that in the 
present model this is ~ because the impurity spin 
ground state is a nonmagnetic singlet. 

IV. SELF ENERGY AND SCATTERING 
MATRIX 

The mean field part of the exchange scattering con- 
tained in Hf^^ is of first order in Jq and therefore a static 
contribution. However the second order dynamic self en- 
ergy due to Hj will exhibit frequency dependence deter- 
mined by the level splittings of the localized spin (Fig[T]). 
The evaluation of the corresponding self energy diagram 
in this figure is described in Appendix |B] The expres- 
sion for the self energy given in Eq. (IBip for general level 
scheme will be simplified to the case which is physically 
relevant here. Firstly we note that the dipolar matrix el- 
ements lead only to two inelastic transitions and only one 
from the ground state. Explicitly (ill^ilO) = M = ^6 
and (±2|S'±| ± 1) = M' = 2. Secondly for T < A the 
M' transition may be neglected since it starts from an 
excited state. Including only the M transition the self 
energy reduces to 

E2(k,r,k'r';zc.„) - iJo/2fp{T) M^^'0^"r") 

k"r" 

^, - l-/(£k"r")+r.B(A) ^ /(6k"r»)+nB(A) l ^^^^ 
iUJn - ek"r" - A ILJn - ek"r" + A J 




FIG. 3. (Color online) Real and imaginary part of self en- 
ergy E2(t^) at T = for ^./E^ = 0.2; /\/Ec = 0.1 and 
(Ao,A) = (0,1) . The singular points given in Eq. H27p at 
high energies close to the surface band edge {H\, H2) and low 
energies close to the Dirac point (L,:, i = 1 — 3) are indicated. 
Full line corresponds to excact result using Eqs. (|23l24p with 
sharp cutoff at Ec and dashed line to numerical integration 
of Eq. (|20p with soft cutoff 7/-Ec = 0.05 and imaginary part 
77 = 0.5 • 10"^ 



where p{T) = 2tanh{A/2T)/[3 - tanh{A/2T)] is the oc- 
cupation difference of singlet ground state and first ex- 
cited doublet and /(w) and nsiij-i) are the Fermi and Bose 
functions respectively. The matrix elements Mj^^, (Ic'V) 
and Matsubara frequencies aj„ are defined in Appendix 
IB] Using Eq. pUj) the total scattering potential for Dirac 
electrons entering into the T- matrix calculation is then 
given by 

f7r'(kk'; iw„) = C^'''(kk') + S2(k, T, kV; ZL^„) (12) 

It is now depending on the energy of conduction elec- 
trons due to the second order exchange self energy. 

The local density of states of conduction electrons at 
the ^ impurity site can be obtained from the expression 

6Nr{uj) = -i^tr[Go(k,c^)r(kk',w)Go(k',w)](13) 

kk' 

where the T-matrix due to the scattering potential Ut 
satisfies the equation 

T(k,k',w) = [/t(k,k';w)-t- (14) 
5];7t(k,k";^)Go(k",c.)T(k",k',c.) 

k" 

The solution of this equation is simplified by the following 
observation. If we neglect the warping term in the Dirac 
Hamiltonian and denote k by its polar coordinates (fc, 0) 
then in Eqs. (|15ll3p Go(fc,a;) and Go{k',uj) depend only 
on the moduli and Uti4'j ^) only on the angles of the 
wave vectors. Then the summation in these equations 
may be done separately for the Green's functions and the 
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FIG. 4. (Color onhne) LDOS for ^ = 0, A = and various 
scattering strengths (Ao,A). Unperturbed D{lj) = -Dol'^l is 
shown as dotted line. Numerical evaluation with r]/Ec — 
0.5 • 10"^ is used. 



T-matrix. Replacing 
the averages 
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-> {A/2n^) J dkd(/) we obtain 
dkGo^k, uj) 



d<j)d<j,'Ut{<j), <j)',uj) 



f{uj) = — / d(^d(^'T{ct),ct)\u) 



(15) 



b' depen- 



Froni Eas. (|AllA2p we can see that in Ut all 
dent terms average to zero. Therefore 

Ut{uj) = [Ua + S2(w)](7o + Uia, 



where (Jq is the unit and az a Pauli matrix. Then Eq.([T5 
may be written as 



(16) 



f{Lo)^Ut{Lo)\l-Ut{u:)Go{^)]-^ 



(17) 



Using helical eigenstates (r = ±1) we may evaluate the 
trace for the LDOS 5N{uj) = J2r SNr{u}). For each he- 
licity channel it is finally given by 



TT 1 - Ul{u})Gl{u;) 



where 



U[{u) = Uo + TUl+t2{u) 



(18) 



(19) 



is the total energy dependent scattering potential at the 
impurity site and {iujn uj + irj) 



2 -JO 

-1- f{er") + nB{A) , f{er")+nB{A) 



(20) 



X 



tUJn - Et 



lUn — Ct 



+ A 



with Et = tE — /i. Here, due to vanishing warping a 
change to integration over the DOS has been used. The 
latter is given by D{E) = DqE with Do = A/{2'kv'^) 
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FIG. 5. (Color online) LDOS at T = for ^ = and 
scattering strengths (Ao,A) = (0,1) for various A. For 
1^1 Ec = 0.2 the thick line is from numerical calculations 
{r}/ Ec — 0.5 ■ 10~^) and the thin line from exact analytical 
results. Unperturbed D{u)) — Do|<^| is shown as dotted line. 
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FIG. 6. LDOS at r = for A/Ec = 0.1, scattering strengths 
(Ao, A) = (0, 1) and various fi. In the lower figure thick line 
represent numerical results {rj/Ec = 0.5 ■ 10~^) and thin line 
the exact analytical results. Arrows indicate position of L2,3 
singularities in self energy due to local spin excitations. 



where A = \s the area of the surface unit cell re- 
lated to the Wigner-Seitz radius fco of the surface BZ by 
fcp — 'in/ A or fco = {2tt / a) / . It is associated with 
the zone-boundary extrapolated cone energy according 
to Eo = vko- The cutoff wave vector is kc = Ec/v <^ fco 
and the cutoff energy Ec <^ Eo is of the order of half 

I 



of the bulk band gap (« 0.25 eV for Bi2Se3) where he- 
lical surface states merge into the bulk states. For nu- 
merical calculations we introduce a soft cutoff function 
/c(a;) = i[l — tanh((i? — Ec)/j)] with 7 denoting the 
width of the cutoff around Ec- This ensures a finite 
width to the bound states around the cutoff energy and 
also improves the convergence of numerical integration. 
Potential scattering and second order self energy may be 
characterized by dimensionless coupling constants Ao and 
A that are, respectively, given by 



Ao = DiEc)Uo 
MT) = -MVoD{Ec)]MT) 



(21) 



Here A depends on temperature due to the occupation 
factor from impurity spin states. Our interest will be 
focused on the case when bound state effects in the LDOS 
appear and these interaction parameters will be chosen 
accordingly in the discussion. The first order exchange 
scattering potential Ui = \ Jo{Sz) will not be considered 
because it vanishes due to the impurity singlet ground 
state, i.e. (OlS'^lO) = and is only weakly induced with 
{Sz) ~ H"^ by a magnetic field. 

The above expressions in Eqs. (I18ll9l20p can now be 
used to caculate the LDOS numerically as function of 
frequency and temperature. A finite imaginary part rj 
for the frequency will be used for the integration. 



V. ANALYTICAL CALCULATION AT T=0 

Instead of presenting immediately numerical results as 
is frequently done within this context, we first develop 
the theory much further to a closed analytical solution 
for T=0 which is possible here due to the restriction to 
isotropic Dirac states with linear quasiparticle DOS. For 
zero temperature we have ^^(A) = 0, /(et) = 1 — 0(/i — 
Er) and p(0)=l. We then may continue to the real axis 
and evaluate the integrals for Gq{uj) = Gq{ijo)' + iGQ{uj)" 
and Yj2{ijj) = S2('^)' + i^2{^)" ■ The former is given by 

Gl {uj)' = Da [-tEc - (c^ + m) In 1 1 - r-^ |] (22) 

UJ + 

Gi{uj)" = -7tDot{lu + ^i)e[T{oJ + ^i)]e[Ec - t{lu + ^)] 

Here 0(x) is the Heaviside function. Likewise we may 
obtain the real and imaginary part of the self energy. 
Defining ^2(0;) = ^M^ J§pi0)^2ioj) - iX/DiEcf)±2iuj) 
we get: 



t2{uJ,fiy = dJ [(c-. - a + Ai)(ln |1 + —^1 - In |1 

1^ LJ — Za 

+ [{uj + A + 



UJ + A' 



-lull 



Ec 



UJ - A- 
Ec 



w + A ■ 



-\)-{u 



A + fi) ln|l + 
A + M)ln|l- 



Er 



OJ + fJ. 
Ec - /i 



OJ — ^ 



"l]e(-M) 

l]e(A*)} (23) 



E2(w, fi)" = -7rDo{(w - A -I- fi)Q{uj - A) [&{uj - A + fi)e{Ec -uj + A-fi)- Q{-uj + A- n)e{Ec +u}- A + fi)] 
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+ (a; + A + /i)e(-cj - A) [6(0; + A + fi)e{Ec - lj - A - ^) - e(-cj 



.)]]{24) 



For chemical potential lying at the Dirac point (/i = 0) 
these expressions simplify to 



+{uj + A) ln|l- 



cj - A' 

Er 



(25) 



and 



±2iuJ, /i)" = -ttDo [{^ - A)e{uj - A)Q{Ec ~ uj + A) 
-{lu + A)e{-uj - A)e(Sc + UJ + A)] (26) 

In this case the self energy fulfils the symmetry relations 
T.2{-(^,fiy = -S2(a;,^)' and S2(-a;,/i)" = I]2(w,/i)". 

The second order self energy has singularities at high 
energies close to the surface band edges uj « ±Ec denoted 
by Hi 2 and at low energies close to the Dirac point de- 
noted by Li_3. The positions of these singularities are 

c.(Li) = (A - M)e(-Ai) - (A + ^^)e{^i) 
uj{L2) = A 
ujiLs) ^ -A 



uj{Hi) = -Ec - (A + fi) 
uj{H2) ^E, + {A-fi) 



(27) 



Because they occur in the denominator of the T- 
matrix they have a direct influence on the positions of 
its bound state poles and therefore on the peak positions 
and anomalies of the LDOS. Only the uj{Li) position de- 
pends explicitly on the chemical potential, it is positive 
for fi < and negative for /i > (A > 0). This singular- 
ity is the weakest because it appears only in the slope of 
the self energy. 



VI. DISCUSSION OF LDOS RESULTS 

In this section wc discuss the results of the theory out- 
lined above. We focus on the effect of second order inelas- 
tic exchange scattering on the LDOS as given by Eq. ([T8|. 
In this expression both real and imaginary part of the 
Green's function and self energy appear. In the whole 
discussion we will ignore the first order exchange scatter- 
ing Ui because it is ^ {Sz) which vanishes for the singlet 
ground state (Fig. [T|) , even in finite field it will increase 
only ^ H^. 

The Green's functions are presented in Fig. [51 In 
Fig.[2ja) the real part for electrons in chiral states t = ±1 
is shown. Full and dashed lines correspond to analyti- 
cal and numerical results respectively. The agreement is 
close keeping in mind that exact results are for a sharp 



I 

cutoff and uj real whereas in the numerical calculation a 
soft cutoff with "f / Ec—Q.05 around the band edge and a 
finite imaginary part i] of the frequency are employed. 
Fig. [2jb) shows the imaginary part which is equal to 
—TrNr{uj)/{DoEc) where Nr is the DOS of the individual 
Dirac cones. 

The T=0 self energy for the general case with both 
fi, A non-zero and different is shown in Fig. [31 Again 
there is an excellent agreement between analytical (full 
lines) and numerical results, even on the details around 
the singular points ^1-3, Hi 2 which may be clearly iden- 
tified. This leads to a confidence that the features in the 
LDOS that originate in the self energy due to inelastic 
scattering will be reliably obtained. We should note that 
the self energy singularities Hi^2 close to the band edges 
and also those at ±Ec in the Green's functions (Fig. ^ 
are to some degree artificial, because they are generated 
by the sharp or soft cutoff at the surface band edges. This 
will also lead to artificial bound states at or outside the 
band edges and corresponding peaks in the model LDOS 
as mentioned in Ref. In reality for w « ±Ec the sur- 
face bands merge into the bulk states which will also show 
a scattering from the impurity potentials. Therefore the 
transition around Ec will be very gradual meaning that 
the high energy singularities will be much weaker. Their 
proper analysis would have to include the surface states 
as well as the bulk states on an equal footing such as de- 
scribed in Refs. [11]. This is beyond the present effective 
surface band model. Therefore in the following when dis- 
cussing the LDOS we will focus on the low energy part 
with uj/Ec < 0.5 which are not much influenced by the 
details of the cutoff procedure and the merging of surface 
and bulk states. 

First we consider the LDOS according to Eq. p^ for 
the symmetric case fi — A = but with variable relative 
strengths of potential (Aq) and exchange (A) scattering. 
The results are shown in Fig. [Hfor various (Aq, A). For 
pure potential scattering Aq ~ 1 > (full line) a strong 
bound state peak develops slightly below the Dirac point 
because the real part 1 — UqGq {oj)' vanishes while the 
imaginary part in Eq. is still small. For the opposite 
sign of the scattering potential Ao ~ — 1 the bound state 
pole develops for the t — + helicity and lies slighly above 
the Dirac point. If we turn on the exchange scatter- 
ing the real part 1 — [/t(a;)Gg"(a;)' vanishes further away 
from the Dirac point because S2(0)' = in the sym- 
metric case. Then the absolute value of the imaginary 
part will be larger. Therefore the bound state signature 
in the LDOS is shifted to lower energy, broadened and 
suppressed as seen in Fig. [H Finally for pure exchange 
scattering A = 1 (Ao — 0) the broad bound state peak 
appears on both sides of the Dirac cone because E2('^)' 
changes sign at a; = 0. In the following we will mostly 
discuss this case because for 3d/4f adatoms the exchange 
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FIG. 7. (Color online) Contour plot of LDOS with Aq = 
0, A(r = 0) = 1. (a) In a; — A plane with chemical potential at 
fi/ Ec — 0.1 ■ Increasing A leads to separation and increasing 
prominence of bound state peaks with asymmetric intensity 
for non-zero /i. (b) In lj — T plane for symmetric case fi — 
and A/Ec = 0.1. 
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FIG. 8. Temperature dependence of LDOS (numerical cal- 
culations for A/Ec = 0.1 with ri/E^ = 0.5 • 10"^). Top: 
fi/Ec = 0. Bottom: /i/E^ = 0.2. 



scattering should dominate due to the effect of strong 
intra-shell correlations. The effect of the inelastic local 
spin excitations on the LDOS for /i = is shown in Fig.[S] 
The Dirac point in the self energy now widens to a finite 
interval [—A, A] according to Eq. (I^Bl) (see also in the 
general case of Fig. [3]). Therefore the bound states are 
pushed symmetrically further away from the Dirac point 
with increasing A/Ec and growing sharper. Remarkably 
however, the self energy singularities L2.3 do no leave a 
signature in the LDOS at ±A. This happens because for 
/i = they are also of the weak type (only in slope) as 
Li is for arbitrary fi. 

The direct signature of inelastic spin excitations finally 
appears in the LDOS if we consider the most general case 
with /X, A different from zero presented in Fig. [6l Then 
ii,2 turn into strong logarithmic singularities (Fig. [3]) 
that can also be seen as anomalies in the LDOS . They 
are stationary in uj as function of fi and show increasing 
prominence for larger /x. This may be understood from 
the terms like (a;-HA+/i)(ln |H-^^|) in Eq. (gS]). For 
a; ~ — A they reduce to /iln | | which is a logarithmic 
singularity at —A whose weight is only finite for non- 
zero fi. Indeed while the effect of this singularity in the 
LDOS is absent in Fig. [5] where /i = its influence grows 
with non-zero /z in Fig. |6l In fact for fi/Ec = 0.2 a 
third broad peak in the spectrum appears right below 
—A which is connected with this inelastic transition. The 
one at A is less prominent but still visible, for ^ < the 
situation would simply be reversed. Since the chemical 
potential may be tuned by applying a gate voltage (2^[25j 
the evolution of these inelastic scattering features in the 
surface LDOS should be observable. The main resonance 
peaks (asymmetric for /.t ^ 0) are already present for 
A = but their distance increases with growing A where 
the asymmetry is also more pronounced as shown in the 
LDOS contour plot in Fig. 

We now discuss the temperature dependence of the ef- 
fects caused by the inelastic local spin excitations on the 
LDOS. There are two temperature effects: Firstly the ef- 



fective exchange coupling A(r) for T > is proportional 
to the occupation factor p(T) which decreases monotoni- 
cally from one at T ^ A to zero for T ^ A. Therefore it 
reduces the size of the self energy due to inelastic scatter- 
ing for larger temperatures. If we assume that potential 
scattering is also present then a change of their relative 
weight as illustrated in Fig. |4] can be achieved by vary- 
ing temperature. Secondly the self energy will depend 
on temperature through the Fermi functions which tend 
to smooth the singularities for increasing temperatures. 
Both effects will suppress the bound state features in the 
LDOS spectrum. The global behaviour is shown in the 
LDOS contour plot of Fig. [TJa for the symmetric case 
fi = 0. While the position of the resonance peaks remain 
almost constant their intensity is rapidly diminished with 
increasing temperature. Cuts at constant temperature 
are shown in Fig.[5]for A/Ec — 0.1. The top panel corre- 
sponds to /i = and T/Ec — 0.1 and the lower panel to 
fi/Ec — 0.2 and T/Ec — 0.05 (full lines). For comparison 
the r = result is also shown (dashed lines). Apparently 
in the case > a moderate temperature is sufficient to 
smooth out the additional bound state structure result- 
ing from the L3- singularity at —A (cf. Fig. [5]). 

We also address the influence of an external magnetic 
fleld on the LDOS spectrum for T = which should be 
considerable a low energies |w| ~ A. In the following we 
consider the easy-plane case {Kq > 0) as before and in 
addition the easy-axis case {Kq < 0) for the local mo- 
ment anisotropy in Eq. ([l} because their LDOS shows 
interesting differences. They correspond to magnetiza- 
tion or susceptibility enhanced parallel or perpendicular 
to the surface respectively. 

First wc discuss the easy-plane case {Kq > 0) as before 
which leads to the level scheme of Fig. [T] We assume a 
symmetric splitting of the | ± 1) excited doublet by the 
fleld into Zeeman levels at A± = AzLS where A = Kq > 
and 6 = qhbH. This corresponds to fleld along z-axis. 
In this case {Sz) is ~ H'^ and therefore the induced Ui 
will be neglected. We also neglect the effect of Landau 
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level splitting on the Dirac electrons for gfisH ^ Ec- 
Therefore only the Zeeman splitting of the impurity spin 
levels enters into the field dependence of the LDOS. It 
is shown in the top part of Fig. [3^ (full line) in com- 
parison to the zero field case (broken line). The main 
effect is a splitting of the zero-field anomaly at —A into 
two well separated features at —A± whose energy differ- 
ence increases linearly with field. In addition the main 
resonance peaks are also slightly modified. 

Now we turn to the complementary easy-axis case 
{Kq < 0) where the level scheme of Fig. [T] will be inverted 
with a doublet | ± 2) ground state and first excited dou- 
blet I ± 1). Both will split linearly in a field along the 
z-axis. Transitions from the ground state to the top- 
most singlet |0) have no dipole matrix element and do 
not appear in the self energy. Furthermore between the 
doublets only one transition | -I- 2) o | -I- 1) is relevant 
at T = 0. Its energy depends on field strength according 
to A+(iJ) = A + S where now A = 3|Xo| > and S as 
before. As a consequence only a single inelastic transi- 
tion appears in the self energy of the easy-axis case. The 
corresponding LDOS is shown in Fig. [9}d and indeed ex- 
hibits only a single additional peak around —A+{H) that 
shifts away from the Dirac point with increasing field. To 
ensure comparison with Fig. IH^ we rescaled the easy-axis 
parameters M' — > M and A — A (or Jq — )• (2/3) Jq and 
lA'ol (l/3)|/^o|). 

The discussion of these two cases shows that the split- 
ting and shifting effect in the LDOS in a magnetic field 
could provide a spectroscopic tool to identify the spin 
states of the adsorbed magnetic impurity even in a very 
dilute limit. In particular the expected different LDOS 
signature for the easy - axis and -plane configurations 
may provide a tool to discriminate between them. 

Finally we estimate the possible influence of the warp- 
ing terms on the inelastic LDOS features. The lowest 
order warping ener gy a ppearing in the massless Dirac 
cone dispersion is [20| e^k — Xkxik^ — 3fc^) . At 
the chemical potential (Fermi wave vector kp) one has 
Cwkp/A ~ (^/A)(^/£'c)^. For parameters used here we 
then have e^^^/A ~ 10^^. Therefore averaging over 
cone states with different warping energies at different 
k-vectors in the LDOS will not wipe out the inelastic 
structures at A discussed above. 

One should note that in reality in Bi2Se3 and Bi2Te3 
the warping of Dirac half-cones is not symmetric and 
furthermore due to self-doping the chemical potential is 
far from the Dirac point 1] where warping effects are 
large. The chemical potential may be shifted close to 
the Dirac point by either carrier doping by excess Se, 
Te or slight chemical substitution of Bi by Ca or Sn [l| 
or application of a gate voltage 0| or a combination of 
both [253. These methods also would have to be applied 
in the present context. 
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FIG. 9. (Color online) Magnetic field effect on LDOS for 
easy-plane and easy-axis anisotropy Ko of impurity spin S 
(a) Ko > : Two transitions to excited impurity states | ± 1) 
with split energies A± = A±5 lead to field induced splitting of 
LDOS inelastic peak (T = 0, fi/E^ = 0.15, A/E^ = 0.1). Thick 
and thin full lines correspond to numerical and analytical re- 
sults respectively. Zero field case is shown as dashed line for 
comparison, (b) Ko < : only one inelastic transition to 
I + 1) excited state is present with shifted energy A+ — A + S 
appearing in the LDOS spectrum. Dashed line as before. 



VII. SUMMARY AND CONCLUSION 

We have investigated the effect of impurity spin states 
split by an anisotropy on the LDOS spectrum of topolog- 
ical insulators. In particular we analyzed a model rele- 
vant for Fe adatoms where the 5 = 2 high spin state leads 
to a nonmagnetic singlet ground state. In this situation 
the usual first order elastic exchange scattering is absent. 
We have demonstrated, however, that in second order of 
the exchange interaction inelastic transitions to excited 
impurity spin states with energy A become important 
and have to be added to the usual nonmagnetic impurity 
potential scattering. The former lead to a dynamic self 
energy contribution in the T-matrix. It has singularities 
associated with the energy scale of inelastic transitions 
that leave its signature in the LDOS firstly by modify- 
ing the position and widths of the bound state features 
in the spectrum already for chemical potential /i = 
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and secondly by introducing new anomalies in the LDOS 
directly connected with the energy of inelastic impurity 
spin excitations. The position of the latter is indepen- 
dent of the chemical potential, however their prominence 
increases with chemical potential moving away from the 
Dirac point. 

Due to thermal occupation the second order self energy 
is rapidly reduced when temperature becomes compara- 
ble to the inelastic transition energy leading to a rapid 
broadening of the bound state features and anomalies as- 
sociated with ±A. We have also suggested that the tem- 
perature dependence of the LDOS spectrum in connec- 
tion with its asymmetry gives direct information about 
the relative weight of pure potential scattering and sec- 
ond order exchange scattering since only the latter de- 
pends on temperature. 

At temperatures T ^ A the inelastic exchange self 
energy anomalies in the LDOS at A may be used as a 
spectroscopic means to determine the surface impurity 
spin splitting microscopically rather than concluding it 
from indirect macrosocpic measurements as described 
in Sec |lTl This could be facilitated by application of a 
magnetic field which leads to a progressive splitting or 
shifting of LDOS anomalies with increasing H, depend- 
ing on the anisotropy type. In this way the inelastic 
exchange processes combined with variation of /i and 
H should make it possible to use LDOS measurements 
as a spectroscopic tool for the spin states of adatoms. 
In particular this would be helpful in a situation when 
several distinct excitation energies A, A'., are accessible 
by exchange scattering which could be realized for 
surface adsorbed magnetic 4f impurities with their 
crystalline electric field splittings. In this case it would 
also be possible to achieve higher excitations energies 
comparable with the full range of the substitution- and 
gate-tuned chemical potential. 

Note added: After completion of this manuscript we 
became aware of a related work on the LDOS of topolog- 
ical insulators in Ref. [2711. 



Appendix A: Matrix elements for helical state 
representation 

In this appendix we give the explicit form of the scat- 
tering matrix elements U'^'^ (kk') and M^Jj (kk') in mean 
field and dynamical exchange Hamiltonian of Eq. pUj) , 
respectively, using the helical basis states. The former 
are given by 



f^o"" (kk) = .g,^^^^ Al 

f/r'(kk') = („,^t^^ ^^'^?+) (Al) 

with A-t = l±expi(0k— 0k') andk'=k+q. The matrix 



elements in the remaining dynamical exchange term "Kj 
are obtained as 



M;j'(kk') = Af-|V;;;(kk') 

F(kk') = P(k')* ® F(k) (A2) 

In general Afisa(25 + l)^x4 matrix where the row 
index a/3 denotes local spin excitations and the column 
index rr' the transition between the helical half-cones. 
The phase factors resulting from helical state represen- 
tation complicate the form of the matrix elements. How- 
ever, as far as the LDOS at the impurity site for Dirac 
cones without warping is concerned, their effect is aver- 
aged out as discussed in Sec. lIVI The explicit form of the 
F-tensor in Eq. (jA2p in terms of the helical states phase 
factors, abbreviating (j) = 4>\^,4>' = 0k' is given by 



V^kk' 



1 








ie"l' 


1 


gi(<t>'+<t>) 


—ie' 




g-i(0'+0) 


1 










1 



icf) 



(A3) 

where the row index is {era') (z-componets of spin) and 
(rr') (helicities) is the column index. 



Appendix B: Evaluation of second order exchange 
self energy 

In this appendix we evaluate the second order self en- 
ergy due to the exchange scattering Hj in Eq. ([5]) . From 
the diagram in Fig. [T]it is obtained as 

S2(k, T, k'r'; iw,,) = {Ja/2fPc ^ M^k' (a/3; k'V") 

a/9;k"T" 

X ^ G{k", r", UJri - ^o) DW, UJm)D{a, UJm + Wo) 

(Bl) 

Here the total matrix element is defined by (see Eq. (jA2p ) 
M^^:{aP;k"T") = M:f{kk") ■ M^y" {k'k") (B2) 

and the Dirac electron and OFF level pseudo-fermion 
Green's functions are given by G{k, r, tj„) = [ia;„— (E'kr— 
/i)]~^ and D{a, ojn) = [itOn — C ~ ^/s]"^ respectively with 
LOn = (2n -|- l)7rT and ujq = 2nQ'!TT denoting fermionic 
and bosonic Matsubara frequencies. Furthermore fi is 
the chemical potential and ^ is a Lagrange parameter to 
project out unphysical localized spin states. The projec- 
tion is performed according to Ref. [2TI with P^[E2(C)] = 
limc^oo^-^eC/^[S2(C)] where Z = Ec^M-E^/T) is 
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the partition function for the spht spin states of Fig. [T] 
We obtain the expression for a general level scheme as 

S2(k,T,kV;zc.„)=(^o/2)2 ^Ck''M;k"T") 

Q:^;k' V 

Pffg [1 - /(£k"T" ) + riB (Agff ] 
ILxJn - (ek"r" + Aq,^) 

Here Ua = Z~^Q:exp(— E'q/T) are the thermal occu- 
pation factors of local spin states and ppa — np — Ua 



their differences. The Fermi and Bose functions are de- 
noted by f{uj) and nB(u}), respectively. Furthermore 
Aq^ = Ea — Ep are their excitation energies and ek"r" = 
EkiiT-"—fJ- quasiparticle energies with respect to the chem- 
ical potential. 
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